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ABSTRACT 

High resolution observations of solar filaments suggest the presence of groups of 
prominence threads, i.e. the fine-structures of prominences, which oscillate coherently 
(in phase). In addition, mass flows along threads have been often observed. Here, we 
investigate the effect of mass flows on the collective fast and slow nonadiabatic mag- 
netoacoustic wave modes supported by systems of prominence threads. Prominence 
fine-structures are modeled as parallel, homogeneous and infinite cylinders embedded 
in a coronal environment. The magnetic field is uniform and parallel to the axis of 
threads. Configurations of identical and nonidentical threads are both explored. We 
apply the T-matrix theory of acoustic scattering to obtain the oscillatory frequency and 
the eigenfunctions of linear magnetosonic disturbances. We find that the existence of 
wave modes with a collective dynamics, i.e. those that produce significant perturba- 
tions in all threads, is only possible when the Doppler-shifted individual frequencies 
of threads are very similar. This can be only achieved for very particular values of the 
plasma physical conditions and flow velocities within threads. 

Subject headings: Sun: oscillations - Sun: magnetic fields - Sun: corona - Sun: 
prominences 



1. Introduction 

Prominences/filaments are fascinating coronal magnetic structures, whose dynamics and prop- 
erties are not well-understood yet. The long life of the so-called quiescent prominences (several 
weeks) suggests that the cool and dense prominence material is maintained against gravity and 
thermally shielded from the much hotter and much rarer solar corona by means of some not well- 
known processes. However, it is believed that the magnetic field must play a crucial role in both the 
support and isolation of prominences. High-resolution observations of solar filaments reveal that 
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they are formed by a myriad of horizontal structures called threads (e.g., iLin et al.l 120051). which 



have been observed in the spines and barbs of both active region and quiescent filaments (|Lin et al 



20081). The width of these fine-structures is typically in the range 0.2 arcsec - 0.6 arcsec, which 
is close to the resolution of present-day telescopes, whereas their lengths are between 5 arcsec 
and 20 arcsec (ILinl 120041) . Threads are assumed to be the basic substructures of filaments and 
to be aligned along magnetic field lines. From the point of view of theoretical modeling, promi- 
nence threads are interpreted as large coronal magnetic flux tubes, with the denser and cooler 
(prominence) region located at magnetic field dips that correspond to the obs erved threads. Al- 
though some theoretical works have attempted to model such structures (e.g- iBallester & Priest 
19891 : ISchmitt & Degenhardtlll995l : iRempel et all 1 19991 : iHeinzel & Anzeiibood) . there are some 
concerns about their formation and stability that have not been resolved yet. 

Small amplitude oscillations, propagating waves and mass flows are some p henomena usually 
observed in prominences and prominence t hreads (see some recent reviews by loiiver & Ballester 
2OO2I : iBallested 120061 : iBanerjee et al.ll2007|) . Periods of small- amplitude prominence oscillations 
cover a wide ran ge from less than a rninute t o several hours, and they are usually attenuated 
in a few periods (|Molowny-Horas et al.l Il999l : iTerradas et al.l 120021) . Focusing on prominence 



threads, some works have detected oscillations and waves in such fine-s tructures (e.g.. 



991I: Yi& Engvold 



199lh and Lin et al 



991I: Linll2004IOkamoto et al.ll2007l : iLin et al.ll2007 [). In particular 



Yi et al. 



Yi et al 



(|2007() suggested the presence of groups of near threads that moved in phase. 



d 

which may be a signature of collecti ve oscillations. On the other hand, mass flows a long magnetic 
field lines have been also detected dzirker et al.lll994[ 1 19981 : lun et all 12001 boosi) . with typical 
flow velocities of less than 30 km i n quiescent prominen ces, although larger values have been 
detected in active region prominences (|Okamoto et al.ll2007|) . Regarding the presence of flows, a 
phenomenon which deserves special attention is t he existence of the so-called coun ter-streaming 
flows, i.e., opposite flows within adjacent threads (|Zirker et al.ll 19981 : iLin et al.ll2003|) . 



Motivated by the observational evidence, some authors have broached the theoretical investi- 
gation of prominence thread oscillations by m eans of the magnetohydrodynamic (MHD) theory in 
the /3 = approximation. First, some works (IJoarder et al.lll997l : iDiaz et al.ll200lL l2003h focused 
on the study of the ideal MHD oscillatory m odes supported by individual nonuniform threads in 
Cartesian geometry. Later, iDfaz et al.l (120021) considered a more representative cylindrical thread 
and obtained more realistic results with respect to the spatial structure of perturbations and the 
behavior of trapped modes. Subsequently, the attention of authors turned to the study of collective 
oscillat ions of groups of threads, and the Cartesian geometry was adopted again for simplicity. 
Hence. iDfaz et al.l (|2005h investigated the collective fast modes of systems of nonidentical threads 
and found that th e only nonleaky mode corresponds to that in which all threads oscillate in spa- 
tial phase. Later, iDfaz & Roberts! (|2006|) considered the limit of a periodic array of threads and 
obtained a similar conclusion. Therefore, these results seem to indicate that all threads within the 



-3- 



prominence should oscillate coherently, even if they have different physical properties. However, 
one must bear in mind that the Cartesian geometry provides quite an unrealistic confinement of 
perturbations, and so systems of more realistic cylindrical threads might not show such a clear 
collective behavior. 

The next obvious step is therefore the investigation of oscillatory r nodes of sys t ems of cylin- 
drical threads. The first approach to a similar problem was done by iLuna et al.l (l2008a|) . who 
considere d a system o f two i dentical, homogeneous cylinders embedded in an unlimited corona. 
Although iLuna et al.l fcOOSah applied their results to coronal loops, they are also applicable to 
prominence threads. These authors numerically found that the system sup ports four trapped kink 



like co llective modes. These results have been analytically re-obtained by Ivan Doorsselaere et al 



J2008I). by conside ring the thin tube approximation and bicylindrical coordinates. Subsequently, 
Luna et al.l (l2008b|) made use of the T-matrix theory of acoustic scattering to study the collective 



oscillations of arbitrary systems of non-identical cylinders. Although the s cattering theory has 
been previously applied in the solar context (e.g.. iBogdan & Zweibellll987l : lKeppenslll994|) . the 
fir st application to the study of normal modes of magnetic coronal structures h as been performe d 
by iLuna et al.l (|2008bh . They concluded that, contrary to the Cartesian case of Diaz et al.l (|2005[ ). 
the collective behavior of the oscillations diminishes when cylinders with nonidentical densities 
are considered, the oscillatory modes behaving in practice like individual modes if cylinders with 
mildly different densities are assumed. 

The present study is based on iLuna et all toOSl h and applies their technique to the investi- 
gation of MHD waves in systems of cylindrical prominence threads. Moreover, we extend their 
model by considering some effects neglected by them. Here, the more general /3 case is consid- 
ered, allowing us to describe both slow and fast magneto acoustic modes. In addi tion, the adiabatic 
assumption is removed and, following previous papers dSoler et al.l 120071 l2008b . the effect of ra- 
diative losses, thermal conduction and plasma heating is taken into account. The detection of mass 
flows in prominences has motivated us to include this effect in our study, and so the presence of 
flows along magneti c field lines is als o considered here. Therefore, the present work extends our 



recent investigation (Soler et al.ll2008 ^. hereafter Paper I), which was focused on individual thread 



oscillations, to the study of collective MHD modes in prominence multi-th read configuration s with 
mass flows. On the other hand, the longitudinal structure of threads (e.g.. iDiaz et al.ll2002 [) is ne- 
glected in the present investigation. For this reason, the effect of including a longitudinal variation 
of the plasma physical conditions within threads should be investigated in a future work. Finally, 
the prominence multi-thread r nodel developed here could be an useful tool for future seismological 
applications (similar to that of iTerradas et al.ll2008b . 



This paper is organized as follows. The description of the model configuration and the math- 
ematical method are given in § |2l Then, the results are presented in § [31 First, the case of two 
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identical prominence threads is investigated in § 13.11 Later, this study is extended to a configura- 
tion of two different threads in § 13.21 Finally, our conclusion is given in § HI 



2. Mathematical method 

Our equilibrium system is made of an arbitrary configuration of homogeneous and unlim- 
ited parallel cylinders, representing prominence threads, embedded in an also homogeneous and 
unbounded coronal medium. Each thread has its own radius, ay, temperature, Tj, and density, py, 
where the subscript j refers to a particular thread. On the other hand, the coronal temperature and 
density are Tc and Pc, respectively. Cylinders are orientated along the z-direction, the xy-plane 
being perpendicular to their axis. The magnetic field is uniform and also orientated along the 
z-direction. By = fiye, being the magnetic field in the j-th thread, and = 5,6, in the coronal 
medium. In addition, steady mass flows are assumed along magnetic field lines, with flow ve- 
locities and directions that can be different within threads and in the corona. Thus, Uy = f/ye, 
represents the mass flow in the 7-th thread, whereas Uc = t/cC, corresponds to the coronal flow. 
For simplicity, in all the following expressions a subscript indicates local equilibrium values, 
while subscripts j or c denote quantities explicitly computed in the j-\h thread or in the corona, 
respectively. 

Such as shown in Paper I, linear nonadiabatic magnetoacoustic perturbations are governed by 
the next equation for the divergence of the velocity perturbation, A = V ■ Vi, 

[rl - [kl + vio) V^] A + AXo^V^A = 0, (1) 
where Tq is the following operator, 

dt "dz 



d d 

To = - + Uo-, (2) 



while v\q = ^ is the Alfven speed squared and A^ is the nonadiabatic sound speed squared. 



a2 = i^o 



(3) 



c^Q = ^ and y being the adiabatic sound speed squared and the adiabatic ratio, respectively. 
Terms with k\\q, tOpo, and ojtq are related to nonadiabatic mechanisms, i.e. radiative losses, thermal 
conduction, and heating (see Paper I for details). Finally, Q.q is the Doppler- shifted frequency 



(iTerra-Homem et al.ll2003l) . 

no = oj- kMo, (4) 
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where co is the oscillatory frequency and is the longitudinal wavenumber. Considering cylindri- 
cal coordinates, namely r, ip, and z for the radial, azimuthal, and longitudinal coordinates, respec- 
tively, we can write A in the following form. 



A = t/r (r, (f) exp {iiot - ik^z) ■ 



(5) 



where the function ij/ (r, ip) contains the full radial and azimuthal dependence. By inserting this last 
expression into equation ([T]), the following Helmholtz equation is obtained. 



where is the Laplacian operator for the r and tp coordinates, and 



(6) 



k^vl^i^-k^A^ 



CjQ - 



(V) 



(8) 



are the radial wave number and the nonadiabatic tube speed squared, respectively. Moreover, 
due to the presence of nonideal terms is a complex quantity. Since nonadiabatic mechanisms 
produce a small correction to the adiabatic wave modes, \''K{ml)\ > \3{ml)\ and the dominant wave 
character depends on the sign of %(ml). Here, we investigate nonleaky modes, which are given by 
5l(m^) < 0. We impose no restriction on the wave character within threads. 



In order to solve equation we consider the technique developed by iLuna et al.l (|200 8b) 
based on the study of wave modes of an arbitrary configuration of cylinders by means of the T- 
matrix theory of acoustic scattering. The novelty with respect to the work of iLuna et al.l (|2008b|) 
is that the method is applied here to solve a Helmhol tz equation for the divergence of the velocity 
perturbation (our eq. [[T]|) whereas iLuna et al.l (l2008b|) considered an equation for the total pressure 
perturbation in the /3 = approximation (their eq. [1]). The present approach allows us to study 
the more general /3 case, therefore slow modes are also described. In addition, nonadiabatic 
effects and mass flows are easily included in our for malism. However, the rest of the technique 
is absolutely equivalent to that of Luna et al.| ( 2008b ). and therefore the reader is refered to their 
w ork for an extensive expl anation of the mathematical technique (see also an equivalent formalism 



in 



Bogdan & Cattaneo 19891) . We next give a brief summary of the method. 



The main difference between our application and that of lLuna et al.l (|2008bl) is in the definition 
of the T-matrix elements. These elements are obtained by imposing appropriate boundary condi- 
tions at the edge of threads, i.e., at |r - Vjl = aj, where is the radial vector corresponding to the 
position of the j-th thread center with respect to the origin of coordinates. In our case, these bound- 
ary conditions are the continuity of the total pressure perturbation, pj, and the Lagrangian radial 
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displacement, = -ivr/^o- Expressions for these quantities as functions of A and its derivative 



are. 



Pt = ipo „3 2 ' ^' 



Expressions for the rest of perturbations are given in Appendix A of Paper I. Thus, in our case the 
T-matrix elements have the following form, 

m,.pj _ fc^v2 J j^^^ (nijaj) J'„ (nicaj) - mjPc (Q^ - ^?v^ ) [nica^ J'„ [mja^ 

T — , (11) 

where Z/^,)^ and J„, are the Hankel function of the first kind and the Bessel function of order m, 
respectively, while the prime denotes the derivative taken with respect to r. Note that the denomi- 
nator of T^^ vanishes at the normal mode frequencies of an individual thread. This can be easily 
checked by comparing it to the dispersion relation of a single thread, see equation (19) of Paper I, 
in which Bessel functions are used instead of Hankel fun ctions. The equiv alence between bot h 



kinds of functions is given in lAbramowitz & StegunI (|l972h . Thus, following iLuna et al.l (|2008bl) . 



the internal i/r (r, (p) field of the j-th thread is, 

DO 

m=-oo 

whereas the external net field is. 



J m=-oo 



where (pj is the azimuthal angle corresponding to the position of the j-th thread center with respect 
to the origin of coordinates, and a^^^ and A/„ are constants. Equations (fT2l) and (fT3l) allow us to 
construct the spatial distribution of A. Subsequently, the rest of perturbations can be obtained. 
Finally, the constants ,„ form a homogeneous system of linear algebraic equations, 

00 

<n + Z Z [^clvj - r,|) = 0, (14) 

ki^j n=-oo 

for -00 < m < 00. Once both integers m and n are truncated to a finite number of terms, the 
non-trivial (i.e., non-zero) solution of system (fT4l) gives us a dispersion relation for the oscillatory 
frequency, a», which is enclosed in the definitions of my and m, . 
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In the next Sections, we apply the method to obtain the oscillatory frequency and the spatial 
distribution of perturbations of the wave modes supported by prominence thread configurations. 
We assume that is real, so a complex frequency is obtained, namely to = to^i + icoi. The imaginary 
part of the frequency appears due to the presence of nonadiabatic mechanisms. The oscillatory 
period, P, and the damping time, td, are related to the frequency as follows, 

P=^, TO = 1. (15) 



3. Results 

3.1. Configuration of two identical threads 

Now, we consider a configuration of two identical threads (see Fig. [T]). Their physical condi- 
tions are typical of prominences (Ti = = 8000 K, pi = p2 = 5 x 10"^' kg m"^) while the coronal 
temperature and density are = 10^ K and pc = 2.5 x 10"'^ kg m"^, respectively. Their radii 
are ai = a2 = a = 30 km, and the distance between centers is J = 4a = 120 km. The magnetic 
field strength is 5 G everywhere. The flow velocity inside the cylinders is denoted by Ui and Ui, 
respectively, whereas the flow velocity in the coronal medium is Uc- Unless otherwise stated, these 
physical conditions are used in all calculations. 



3.1.1. Wave modes in the absence of flow 



First, we consider no flow in the equilibrium, i.e. Ui = U2 = Uc = 0. We fix the longitudinal 
wavenumber to ha = 10~^, which corresponds to a wavelength within the typically observed 
range. In addition to the four kink modes described by Luna et al. ( 2008a[) . i.e. the 5 v, A^-, Sy, and 
Ay modes, where 5 or A denote symmetry or antisymmetry of the total pressure perturbation with 
respect to the jz-plane, and the subscripts refer to the main direction of polarization of motions, we 
also find two more fundamental collective wave modes (one symmetric and one antisymmetric) 
m ainly polari z ed alon g the z-direction, which we call S ^ and A^ modes following the notation 
of Luna et al. ( 2008al). These new solutions correspond to slow modes which are a bsent in the 



investig ation of iLuna et al.l (|2008ah due to their /3 = approximation. As was stated by iLuna et al. 



(|2008b|) . a collective wave mode is the result of a coupling between individual modes. So the reader 
must be aware that in the present work we indistinctly use both expressions, i.e., collective modes 
and coupled modes, to refer to wave solutions whose perturbations have significant amplitudes in 
both threads. 



The total pressure perturbation field, pj, and the transverse Lagrangian displacement vector- 
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field, corresponding to the six fundamental modes are displayed in Figure[2l On the other hand, 
Figure [3] displays a cut of the Cartesian components of the Lagrangian displacement {^x, ^y, and 
at J = 0, again for these six solutions. For simplicity, only the real part of these quantities are 
plotted in both Figures, since their imaginary parts are equivalent. One can see in Figure [2] that 
the amplitude of the longitudinal (magnetic field aligned) Lagrangian displacement, of the S ^ 
and Aj. modes is much larger than the amplitude of transverse displacements, and ^y, such as 
corresponds to slow modes in/3 < 1 homogeneous media, while the contrary occurs for the S x, Aj^, 
Sy, and Ay fast kink solutions. 

Next, Figure displays the ratio of the real part of the frequency of the four kink solutions 
to the frequency of the individual kink mode, tOk (from Paper I), as a func tion of the distance 
between the center of cylinders, d. This Figure is equivalent to Figure 3 of iLuna et al. I (|2008a) 
and, in agreement with them, one can see that the smaller the distance between centers, the larger 
the interaction between threads and so the larger the separation between frequencies. On the other 
hand. Figure |4j5 shows the ratio of the damping time to the period of the four kink modes as a 
function of d. We see that the damping times are between 4 and 7 orders of magnitude larger 
than their corresponding periods. Therefore, dissipation by non-adiabatic mechanisms cannot be 
responsible for the o bserved damping tirn es of transverse thread oscillations, as was pointed out 
in Paper I. Recently, lArregui et al.l (|2008l) found that the mechanism of resonant absorption can 
provide kink mode damping times compatible with those observed. 

Regarding slow modes. Figure [5^ displays the ratio of the real part of the frequency of the S ^ 
and A^ solutions to the frequency of the individual slow mode, oj., (from Paper I). One can see that 
the frequencies of the and A, modes are almost identical to the individual slow mode frequency, 
and so the strength of the interaction is almost independent of the distance between cylinders. This 
is consistent with the fact that transverse motions (responsible for the interaction between threads) 
are not significant for slow-like modes in comparison with their longitudinal motions. Therefore, 
the S z and A^ modes essentially behave as individual slow modes, contrary to kink modes, which 
display a more significant collective behavior. Finally, Figure [5h shows td/P corresponding to the 
Sy and Ay solutions versus d. One sees that both slow modes are efficiently attenuat ed by non 



adiabaticmechanisms, with td/P x 5, which is in agreement with previous studies dSoler et al 



20071 12008 ) and consistent with observations. 



3.1.2. Effect of steady mass flows on the collective behavior of wave modes 

The aim of the present section is to assess the effect of flows on the behavior of collective 
modes. With no loss of generality, we assume no flow in the corona, i.e. Uc = 0. On the other 
hand, the flow velocities in both cylinders, namely Ui and U2, are free parameters. We vary these 
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flow velocities between -30 km s"^ and 30 km s~\ which corr espond to the range of typically 
observed flow velocities in filament threads (e.g. jLin et al.ll2003|) . These flow velocities are below 
the critical value that determines the apparition of the Kelvin-Helmholtz instability (see details in 
Holzwarth et al. In our configuration, a positive flow velocity means that the mass is flowing 

towards the positive z-direction, wh ereas the contrary is for negative flow velocities. From Paper I 
(see also lTerra-Homem et al. I I2OO3I) we know that the symmetry between waves whose propagation 
is parallel (cur > 0) or anti-parallel {co^ < 0) with respect to magnetic field lines is broken by the 
presence of flows. Hence, we must take into account the direction of wave propagation in order 
to perform a correct description of the wave behavior. Following Paper I, we call parallel waves 
those solutions with cor > 0, while anti-parallel waves are solutions with cor < 0. 

We begin this investigation with transverse modes. First, we assume Ui = 20 km s"' and 
study the behavior of the oscillatory frequency when U2 varies (see Fig. (6]). Since frequencies 
are almost degenerate and, therefo r e, alrr iost indiscernible if they are plotted together, we use the 
notation of Ivan Doorsselaere et al.l (|2008|) and call low-frequency modes the S x and Ay solutions, 
while high-frequency modes refer to Av and Sy solutions. In addition, we restrict ourselves to 
parallel propagation because the argumentation can be easily extended to anti-parallel waves. To 
understand the asymptotic behavior of frequencies in Figure [6l we define the following Doppler- 
shifted individual kink frequencies: 



^kl = OJk + U2K. 

Since Ui is fixed, Q.k\ is a horizontal line in Figured whereas 0.k2 is linear with Uj- 



(16) 
(17) 



Three interesting situations have been pointed by means of small letters from a to c in Fig- 
ure [6l Each of these letters also corresponds to a panel of Figure |7] in which the total pressure 
perturbation field of the S x mode is plotted. The three difl'erent situations are commented in detail 
next (remember that in all cases Ui =20 km s"^): 



• a) U2 = -10 km s"^ {U2 < Ui). This corresponds to a situation of counter-streaming flows. 
From Figured we see that the frequency of low-frequency modes is close to Qti, whereas 
that of high-frequency solutions is near Qki ■ Thus, these solutions do not interact with each 
other and low-frequency (high-frequency) solutions are related to individual oscillations of 
the second (first) thread. This is verified by looking at the total pressure perturbation field in 
Figure |7^, that shows that only the second thread is significantly perturbed. Therefore, for 
an external observer this situation corresponds in practice to an individual thread oscillation. 

• b) U2 = 20 km s"^ {U2 = Ui). The flow velocities and their directions are equal in both 
threads. In such a situation, low- and high-frequency modes couple. At the coupling, an 
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avoided crossing of the solid and dashed lines is seen in Figure [6l Because of this cou- 
pling solutions are related no more to oscillations of an individual thread but they are now 
collective and, for this reason, Figure |7}) shows a significant pressure perturbation in both 
threads. 

• c) U2 = 27 km s"' {U2 > U\). This case is the opposite to situation a). Therefore, in practice 
the present situation corresponds again to an individual thread oscillation. 

Our argumentation is supported by Figure [H which displays the amplitude of the transverse 
Lagrangian displacement, at the center of the second thread as a function ofU2, for parallel and 
anti-parallel kinklike waves. The displacement amplitude at the center of the first thread is always 
set equal to unity. The three previously commented situations have been pointed again in Figure [81 
We clearly see that the displacement amplitude is only comparable in both threads, and so their 
dynamics is collective, when their flow velocities are similar. 

Next we turn our attention to slow modes. The behavior of the S ^ and modes is like that of 
low- and high-frequency kinklike solutions, so we comment them in short for the sake of simplicity. 
Sy and Ay solutions can only be considered collective when the flow velocity is the same in both 
threads because, in such a case, the S ^ and A, modes couple. If difl^erent flows within the threads 
are considered, the S ^ and A, slow modes lose their collective aspect and their frequencies are very 
close to the Doppler-shifted individual slow frequencies, 

ai = CO, + Uik„ (18) 

^2 = CO, + U2K. (19) 

Then, the S ^ and A^ solutions behave like individual slow modes. It is worth to mention that the 
coupling between slow modes is much more sensible to the flow velocities in comparison with fast 
modes, and the S ^ and Ay solutions quickly decouple if Ui and U2 slightly difi'er. An example of 
this behavior is seen in Figure[9l which displays the total pressure perturbation field of the Ay mode 
for Ui - U2 = 10"^ km s~'. Although the difference of the flow velocities is insignificant, one can 
see that in this situation the A, mode essentially behaves as the individual slow mode of the second 
thread. 

The main idea behind these results is that fast or slow wave modes with a collective appear- 
ance (i.e., modes with a similar displacement amplitude within all threads) are only possible when 
the Doppler-shifted individual kink (eqs. f\M and IfTTll ) or slow (eqs. [[TSlI and |[T9l ) frequencies are 
similar in both threads. In a system of identical threads, this can only be achieved by considering 
the same flow velocities within all threads, since all of them have the same individual kink and 
slow frequencies. However, if threads with difl^erent physical properties are considered (i.e., differ- 
ent individual frequencies), it is possible that the coupling may occur for different flow velocities. 
This is explored in the next section. 
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3.2. Configuration of two threads with different physical properties 

Now we consider a system of two nonidentical threads and focus first on kink modes. From 
§ 13.11 we expect that collective kink motions occur when the Doppler-shifted individual kink fre- 
quencies of both threads coincide. The relation between flow velocities U\ and U2 for which the 
coupling takes place can be easily estimated from the equation, 

(x>k\ + Uih ~ (x)k2 + U2K. (20) 

Note that the individual kink frequency of each thread is diff'erent, i.e., oj^i + o^ki- Then, the 
relation between flow velocities at the coupling is, 

U,-U,^ <^^^^Z^_ (21) 

K 

This last expression can be simplified by considering the approximate expression for the kink 
frequency in the long-wavelength limit. 



<X)ki ~ + J- ]—VAiK « ± ^VAiK, (22) 

1 + pjpi 

for i = 1,2, where the -i- sign is for parallel waves and the - sign is for anti-parallel propagation. 
Then, one finally obtains, 

t/l-[/2«±V2(VA2-VAl), (23) 

where the meaning of the -I- and - signs is the same as before. In the case of identical threads, 
vai = va2 and so t/i - t/2 = 0. Thus the flow velocity must be the same in both threads to obtain 
collective motions, as we concluded in § 13.11 An equivalent analysis can be performed for slow 
modes and one obtains, 

f/i-[/2«±(c,2-c,i). (24) 

In general, the coupling between slow modes occur at different flow velocities than the coupling 
between kink modes. This makes difficult the simultaneous existence of collective slow and kink 
solutions in systems of non-identical threads. 

Next, we assume a particular configuration of two non-identical threads in order to verify 
the later argumentation. Thread radii are ai = 30 km and ^2 = 45 km, whereas their physical 
properties are = 8000 K, pi = 5 x 10"" kg m'^, and T2 = 12000 K, p2 = 3.33 x 10"^' kg m ^ 
Coronal conditions are those taken in § 13.11 (i.e.. Tc = 10^ K, pc = 2.5 x 10 kg m"^). The 
magnetic field strength is 5 G everywhere and the distance between the thread centers is d = 
120 km. We assume t/i = 10 km s"^ In the present case, also four kinklike solutions are present, 
which are grouped in two almost degenerate couples. For this reason, we again refer to them 
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as low- and high-frequency kink solutions. Their frequency as a function of U2 is displayed in 
Figure [lOl At first sight, we see that solutions couple for a particular value of U2, as expected. 
Applying equation (|23] ). and considering that Vai = 63.08 km s~' and v^i = 77.29 km s~\ we 
obtain U\ - U2 ~ +20.10kms^\ and since Ui = 10 km s"\ we get U2 ~ -10.10 km s~' for 
parallel propagation. We see that the approximate value of t/2 obtained from equation (|23]) is in 
good agreement with Figure \T0\ In addition, we obtain that for parallel propagation collective 
dynamics appear in a situation of counter- streaming (opposite) flows. This result is of special 
relevance becau se counter-streaming flows have been detected in prominences dZirker et al.lll998 : 



Lin et al.l 120031) and might play a crucial role in the collective behavior of oscillations. On the 



contrary, in the anti-parallel propagation case we obtain U2 « 30.10 km s"' from equation (|23]) . 
and so both flows are in the same direction and large a quite value of U2 is obtained in comparison 
with the parallel propagation case. 

Regarding slow modes, considering that c.vi = 11.76 km s~^ and Cs2 = 14.40 km s"\ equa- 
tion (|24l) gives U2 « 7.36 km s"' for parallel propagation and U2 « 12.64 km s"^ for anti-parallel 
waves. Note that in our particular example the flow velocities needed for the coupling situation are 
realistic and within the range of typically observed velocities. However, if threads with very differ- 
ent physical properties and, therefore, with very diff"erent Alfven and sound speeds are considered, 
the coupling flow velocities could be larger than the observed values. This means that the condi- 
tions necessary for collective oscillations of systems of threads with very different temperatures 
and/or densities may not be realistic in the context of solar prominences. 



4. Conclusion 



In this work, we have assessed the effect of mass flows on the collective behavior of slow 
and fast kink magnetosonic wave modes in systems of prominence threads. We have seen that the 
relation between the individual Alfven (sound) speed of threads is the relevant parameter which 
determines whether the behavior of kink (slow) modes is collective or individual. In the absence 
of flows and when the Alfven speeds of threads are similar, kink modes are of collective type. 
On the contrary, perturbations are confined within an individual thread if the Alfven speeds differ. 
In the case of slow modes, the conclusion is equivalent but replacing the Alfven speeds by the 
sound speeds of threads. On the other hand, when flows are present in the equilibrium, one can 
find again collective motions even in systems of non-identical threads by considering appropriate 
flow velocities. These velocities are within the observed values if threads with not too different 
temperatures and densities are assumed. However, since the flow velocities required for collective 
oscillations must take very particular values, such a special situation may rarely occur in real 
prominences. 
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Therefore , if coherent oscillations of groups of threads are observed in prominences (e.g., 



Lin et al.ll2007() . we conclude that either the physical properties and flow velocities of all oscillat- 
ing threads are quite similar or, if they have diff'erent properties, the flow velocities within threads 
are the appropriate ones to allow collective motions. From our point of view, the first option is the 
most probable one since the flow velocities required in the second case correspond to a very pe- 
culiar situation. This conclusion has important repercussions for future prominence seismological 
applications, in the sense that if collective oscillations are observed in large areas of a prominence, 
threads in such regions should possess very similar temperatures, densities, and magnetic field 
strengths 

Here, we have only considered two-thread systems, but the method can be applied to an 
arbitrary multi-thread configuration. So, the model developed here could be used to perform seis- 
mological studies of large ensembles of prominence threads if future observations provide with 
positions and physical parameters of such systems. 
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Fig. 1. — Scheme in the xj-plane of the model considered in § 13.11 The z-axis is perpendicular to 
the plane of the figure and points towards the reader. 
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Fig. 2. — Real part of the total pressure perturbation field (contour plot in arbitrary units) and the 
transverse Lagrangian displacement vector-field (arrows) plotted in the xj-plane corresponding to 
the wave modes a) S b) A^, c) Sy, d) Ay, e) S ^, and /) A^ in the absence of flows for a separation 
between threads d = 4a and a longitudinal wavenumber k^a = 10"^. The location of prominence 
threads is denoted by dotted circles. 
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Fig. 3. — Cut at y = 0, z = of the real parts (in arbitrary units) of the Cartesian components of 
the Lagrangian displacement: (solid line), (dotted line), and ^, (dashed line), corresponding 
to the wave modes a) S x, b) A^, c) Sy, d) Ay, e) S and /) A, for the same conditions of Figure[2l 
The shaded regions show the location of threads. Note that neither ^y nor are continuous at the 
edges of threads. 
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Fig. 4. — a) Ratio of the real part of frequency, wr, of the S x (solid line), (dotted line), Sy 
(triangles), and Ay (diamonds) wave modes to the frequency of the individual kink mode, cj^, as 
a function of the distance between centers, b) Ratio of the damping time to the period versus the 
distance between centers. Linestyles are the same as in panel a). 
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Fig. 5. — a) Ratio of the real part of frequency, a>R, of the S ^ (solid line) and (dotted line) wave 
modes to the frequency of the individual slow mode, oj^, as a function of the distance between 
centers, b) Ratio of the damping time to the period versus the distance between centers. Linestyles 
are the same as in panel a). 
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Fig. 6. — Ratio of the real part of the frequency, a»R, to the individual kink frequency, to^, as a 
function of U2 for U\ = 20 km s"^ The solid line corresponds to parallel low-frequency modes 
{S X and Ay) while the dashed line corresponds to parallel high-frequency solutions (A^ and Sy). 
Dotted lines correspond to the Doppler- shifted individual kink frequencies of the threads, 0.^1 and 
0.k2- The small letters next to the solid line refer to particular situations studied in the text. 
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Fig. 7. — Real part of the total pressure perturbation field (contour plot in arbitrary units) and the 
transverse Lagrangian displacement vector-field (arrows) plotted in the xy-plane corresponding to 
the parallel S x mode for a) U2 = -10 km s"\ b) U2 = 20 km s^\ and c) U2 = 27 km s"^. In all 
cases Ui =20 km s~K 
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Fig. 8. — Amplitude (in arbitrary units) of the transverse Lagrangian displacement, at the center 
of the second thread as a function of Uj. The amplitude displacement in the first thread center is 
fixed to unity (dotted line) and the flow velocity is t/i = 20 km s~'. Solid lines correspond to low- 
frequency kink modes (5 ^ and Ay), while dashed lines correspond to high-frequency kink solutions 
(A v and Sy). Lines without simbols are for parallel waves whereas diamonds indicate anti-parallel 
propagation. Small letters refer to the particular situations discussed in the text. 
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Fig. 9. — Real part of the total pressure perturbation field plotted in the xy-plane corresponding to 
the parallel A J. slow mode for a difference between flow velocities of t/i - t/2 = 10"^ km s"'. 
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Fig. 10. — Same as Figure[6]but for the particular configuration studied in § 13.21 



